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Abstract 

The Aharonov-Bohm effect is considered as a scattering event with 
nonrelativistic charged particles of the wavelength which is less than the 
transverse size of an impenetrable magnetic vortex. The quasiclassical 
WKB method is shown to be efficient in solving this scattering problem. 
We find that the scattering cross section consists of two terms, one describ- 
ing the classical phenomenon of elastic reflection and another one describ- 
ing the quantum phenomenon of diffraction; the Aharonov-Bohm effect is 
manifested as a fringe shift in the diffraction pattern. Both the classical 
and the quantum phenomena are independent of the choice of a boundary 
condition at the vortex edge, providing that probability is conserved. We 
show that a propagation of charged particles can be controlled by altering 
the flux of a magnetic vortex placed on their way. 
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1 Introduction 



A comprehensive description of electromagnetism in classical theory is given in 
terms of the electromagnetic field strength acting locally and directly on charged 
matter. This is not the case in quantum theory where a state of charged mat- 
ter is influenced by electromagnetic field even in the situation when the region 
of the nonvanishing field strength does not overlap with the region accessible 
to charged matter; the indispensable condition is that the latter region be non- 
simply-connected. In particular, a magnetic field of an infinitely long current- 
carrying solenoid is made impenetrable to charged matter; such a field configura- 
tion may be denoted as an impenetrable magnetic vortex. Perhaps, W. Ehrenberg 
and R. E. Siday were the first to discover theoretically the purely quantum effect 
which is due to the influence of the vector potential of the impenetrable magnetic 
vortex on the motion of a charged particle, and they proposed an interference 
experiment to detect this pp. Ten years later the effect was rediscovered by 
Y. Aharonov and D. Bohm who made an important addition by predicting also 
an effect which is due to the influence of the electric scalar potential [2j. Immedi- 
ately, in a year, the first interference experiment to verify the magnetic effect was 
performed [3] , and this caused a burst of interest opening a long way of an exten- 
sive research both in theory and in experiment, see, e.g., [U El El CO El El HOI ED] 
and references therein; the effect conventionally bears the name of the authors of 
[2] . Although a great progress in understanding the Aharonov-Bohm ( AB) effect 
has been achieved since the time of foundational works pfl El E] , one may agree 
with the authors of [12J that "the investigation and exploitation of the AB effect 
remain far from finished" . 

It should be noted that the concern of paper [I] was in the elaboration of 
the consistent theory of refraction of electron rays in the framework of electron 
wave optics. Following thoroughly their flow of reasoning, the authors arrive at 
the conclusion that the medium in a spatial region outside a magnetic vortex is 
anisotropic and that there exist wave-optical phenomena which are due to the 
influence of the vector potential of the vortex. Finally, in the very end of the 
paper, they suggest an interference experiment to verify their findings. 

An approach of paper [2] was quite different: the authors intuitively, with 
the use of few rather general and obvious formulas (which are even left unnum- 
bered in the text), arrive at their suggestion of interference experiments already 
in the beginning of the paper. A powerful theoretical background behind this 
heuristic suggestion is hereafter given in Section 4. "Exact solution for scattering 
problems" (which exceeds in its volume two preceding sections). The background 
relies directly on the formalism of quantum mechanics, namely on the Schrodinger 
wave equation, which allows the authors to consider their effect as a scattering 
event. Perhaps, this was a reason, why paper [2J, unlike paper [TJ, had received 
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an immediate response from the physics community. 

In this respect it looks rather paradoxical that, whereas interference exper- 
iments aiming at the verification of the AB effect were numerous (see [H]), a 
scattering experiment aiming at the same was never performed. One of the mo- 
tivations of the present paper was to provide an explanation for this paradox. 

Let us commence by recalling the famous formula obtained by the authors 
of [2] for the scattering differential cross section per unit length of the magnetic 
vortex (see (21) and (22) in 0): 

dcr( AB ) h sin 2 [e$(2^c)- 1 ] 

dz&v ~ 2^) sin%/2) ' ^ 

where p is the momentum of a scattered particle, e is its charge, p is the scattering 
angle (with (p = corresponding to the forward direction) and $ is the flux of 
the magnetic vortex directed along the z-axis. A characteristic feature of (1) is a 
periodic dependence in the magnetic flux with the period equal to the London flux 
quantum, 2nhc/\e\. Also, one could suggest that (1) describes a purely quantum 
effect, since it vanishes in the formal limit of h — > 0, or, more physically, in 
the limit of large enough values of the particle momentum, p — > oo. However, 
result (1) corresponds to the idealized case of the magnetic vortex of the vanishing 
transverse size. Since any real magnetic vortex (either that produced by a current- 
carrying solenoid, or that inside a magnet) is of a nonzero transverse size, the 
quotient of squared sines in (1) in the real case becomes a certain function of 
dimensionless variable pr c /h, where r c is a length characterizing the transverse 
size of the vortex. The transverse size effects were taken into account in [TJ, and 
results were found to converge with (1) in the case of pr c hr x — > 0. As to the case 
of pr c h~ x ^> 1, it was only scarcely and insufficiently considered, see [HI E], and, 
therefore, let us dwell on this case at more length. 

It should be noted in the first place that the smallest possible values of r c 
can be of order 1 fim (ferromagnetic whiskers enclosed in screening shells) [llj, 
while, in a scattering event, the largest possible wavelengths, A = 2ttH/p, of the 
lightest charged particle, electron, are of order 0.1 nm (slowly moving electrons 
with energies 10-100 eV, e.g., such as those in famous diffraction experiments 
by C. J. Davisson and L. H. Germer in 1927, Nobel prize in 1937). Hence, 
quantity pr c h~ x takes values of order 10 4 — 10 5 at least, and formula (1) which 
might be supposed to approximate the case of pr c hr x <C 1 has no relation to 
scattering of real particles. Albeit this formula looks essentially quantum, the 
"quantum" factor of h/p appears just due to dimensional reasons, because r c = 0; 
it is a somewhat more explicating to write h/p as product r c x (prc/h) -1 , where 
the second (dimensionless) factor is vanishingly small for all charged particles 
existing in nature. This was certainly understood by the authors of [2], who 
never proposed cross section (1) for the experimental verification. 
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It may seem quite reasonable to apply the quasiclassical (called also semiclas- 
sical) method of G. Wentzel, H. A. Kramers and L. Brillouin (WKB) [T3l[TllT5] 
to solving the problem of scattering of real charged particles by an impenetrable 
magnetic vortex: the size of the interaction region may be estimated as being of 
order of r c and larger, then condition r c A _1 ^> 1 allows one to hope that the use 
of the WKB method will be adequate, see, e.g, [161 HZ] - However, a scepticism 
immediately arises: how can a purely and essentially quantum effect as is the AB 
one be deduced with the help of a quasiclassical (semiclassical) method which 
is commonly believed to give nothing else but the yield of classical theory? In 
this respect it should be noted that the findings of the authors of [6j [8] indicate 
that scattering off an impenetrable magnetic vortex in the case of prjir 1 ^> 1 is 
reduced to classical scattering of point particles off an impenetrable tube with no 
flux. 

Nevertheless, in spite of the scepticism substantiated by the results of [HI E], 
we shall proceed by applying the WKB method to the solution of the scattering 
AB problem in the present paper. Let us recall that interference experiments 
aiming at the verification of the AB effect were performed with rather energetic 
quasiclassical electrons. We are considering electrons of the same (as in the 
interference experiments) and even smaller (down to 10-20 eV) energies and study 
direct scattering of such particles on an impenetrable magnetic vortex, using the 
WKB method. Whereas the authors of j6j [7J [8] employed the Dirichlet boundary 
condition to ensure the vortex impenetrability, we go further by employing the 
most general (unorthodox in the language of [6]) boundary condition that is 
required by the basic concepts of probability conservation and self-adjointness of 
the Hamilton operator. In the course of our study we shall find, in addition to 
the classical effect of [HI E] , a purely quantum effect that can be denoted as the 
scattering AB effect with real particles. 

In the next section we consider a quantum-mechanical charged particle in the 
background of an impenetrable magnetic vortex and choose a boundary condi- 
tion for the particle wave function at the vortex edge. The WKB method is 
used to obtain the partial waves of the solution to the Schrodinger equation in 
Section 3. The scattering amplitude and cross section are obtained in Section 4. 
The conclusions are drawn and discussed in Section 5. We relegate some details 
of calculation of the scattering amplitude to Appendix A, while the use of the 
Hamilton- Jacobi equation of motion to derive the classical effect is outlined in 
Appendix B. 
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2 Self-adjointness of the Hamilton operator and 
the Robin boundary condition 

Defining a scalar product as ( X , x) — J d 3 xx*X, we get , using integration by 
parts, 



1 / e \ 2 

-— (hd-i-A) +v 

2m V c / 



an 



ftc? — i-A ) x 



(hd + i e -A) r] x) = (H% x) ~ rtj da- ■ u, (2) 



i-m 

where dfl is a twodimensional surface bounding the threedimensional spatial 
region Q, 

H = H, = -^( hd - i l A ) 2 + v < 3 > 

is the formal expression for a general quantum-mechanical Hamilton operator 
in an external electromagnetic field (electric scalar potential is included in the 
general potential energy V), and 



u 



i r [( hd ~ ' l ~c A ) A -[{ hd+ { - c A ) *i • 

Operator H is Hermitian (symmetric), 



if 



dcr ■ u = 0. 



(4) 



(5) 



(6) 



The latter condition can be satisfied in various ways by imposing different bound- 
ary conditions for X an d X- However, among the whole variety, there may exist 
a possibility that a boundary condition for x is the same as that for x; then the 
domain of definition of W (set of functions x) coincides with that of H (set of 
functions x)i an d operator H is called self-adjoint. Whether such a possibility 
exists is in general determined by the Weyl - von Neumann theory of self-adjoint 
operators, see, e.g., [18]. The action of a self-adjoint operator results in func- 
tions belonging to its domain of definition only, and, therefore, a multiple action 
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and functions of such an operator (for instance, the evolution operator) can be 
consistently defined. 

In the present case, the problem of self-adjointness of operator H (3) is re- 
solved by imposing a boundary condition in the form 



sin(p7r) n 
sin(p7r) n • 



le 



d - —A ) x + cos(pvr) 



d 



le 



n • x 

1 



-X 



-Ajx + cos( p7 r) n x 



X 



x€<9fi 

= o, 



(7) 



xedQ 



where n is the internal unit normal to boundary dfl and p is the real parameter 
(self-adjoint extension parameter). Defining current 



J = -d;k[( M - i 3 A )*H('" ,+ 1 A ) **]*}• 



(8) 



we note that self-adjointness condition (7) results in the vanishing of the normal 
component of the current at the boundary: 



n ■ j|xe<9n = 0. 
If X is a solution to the Schrodinger equation, 



(9) 



(10) 



then current (8) is the probability current density satisfying the continuity equa- 
tion 

!rlx 2 | + 0'j = O. 



dV 



(11) 



Thus self-adjointness of the Hamilton operator ensures probability conservation, 

|/^lxl 2 = o. 

n 

In the following we consider a quantum-mechanical charged particle in the 
background of a static magnetic field, hence V = 0. The magnetic field is confined 
to an infinite tube which is classically impenetrable to the charged particle; thus, 
spatial region Q is an infinite threedimensional space with the exclusion of the 
tube containing the magnetic flux lines, and boundary surface dQ is a surface of 
the tube. Assuming the cylindrical symmetry of the flux tube, we use cylindrical 
coordinates x = (r, ip, z) with the z-axis coinciding with the axis of the tube and 
choose the components of the vector potential in Q in the form 



0, 



(12) 
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where $ is the total magnetic flux confined in the tube. Then the Hamilton 
operator takes the form 



H = — 



2m 



-d r rd r + \ 



d in - 



ie$ V 
2nhc) 



(13) 



and the solution to Schrodinger equation (10) with H (13) is 
X(t, x) = exp(—iEth~ + Yp z zfr x )^{r, if), 
where ip(r, if) is the solution to equation 



r r r r 2 \ ^ 2nhc J K 2 



ip(r,if) = 0, 



and obeys the boundary condition, see (7), 



[sm.(fm)rd r 'il){r ) <f) + cos(/97r)^(r, if)] | r= 



0; 



(14) 



(15) 



(16) 



r c is the tube radius, E = (p^+p 2 )(2m) _1 is the energy of the stationary scattering 
state, p z and p are the momenta in the longitudinal and the radial transverse 
directions, respectively. One could recognize that (16) is known as the Robin 
boundary condition; the case of p = corresponds to the Dirichlet condition 
(perfect reflectivity of the boundary) and the case of p = 1/2 corresponds to the 
Neumann condition (absolute rigidity of the boundary). It should be emphasized 
that parameter p is in general dependent on z and if. Thus, the "number" of 
self-adjoint extension parameters is infinite, moreover, it is not countable but is 
of power of a continuum. This distinguishes the case of an extended boundary 
from the case of an excluded point (contact interaction) when the number of self- 
adjoint extension parameters is finite, being equal to n 2 for the deficiency index 
equal to {n,n} (see, e.g., [T9]). 

The motion of the particle in the longitudinal direction is free, and our aim is 
to determine the particle motion in the orthogonal plane, i.e. to find a solution 
to (15) and (16). 



3 WKB method 

Let us start with the decomposition of wave function ip{r, if) into partial waves 

ij(r,f) = J2e lnip a n R n (r), (17) 
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where Z is the set of integer numbers. The radial component of each partial wave 
satisfies equation 



[A r + P*(r)fr- 2 ]Rn(r) = 0, 



where A r = r l d r rd r and 



P n {r) = \J p 2 — h 2 (n — fi) 2 r~ 



(19) 



fi = e$(27r7ic)~ 1 . Note that ip(r,(p) (17) is periodic in the value of flux $ with 
the period equal to the London flux quantum. 

We present the radial component of the partial wave as R n {r) = exp [^E n (r 
where £ n (r) satisfies the Ricatti-type equation 



^A r S n -l(9 r S n ) 2 + -ip n 2 = 0. 



(20) 



Expanding £ n (r) into series 



1=0 



h 



ipr f 



(21) 



and collecting terms of the same order in H, we obtain the system of equations: 

2 



<9 r En 



(0) 



<9 r £ n 



(0) 



<9 r £ n 



(0) 



<9 r £ n 



(i) 



<9 r £ n 



(2) 



pr c A r £ n ^ , 



pr c A r £ n 



(i) 



d r Yj n 



(i) 



(0) 



d r Yj n 



(0) 



(21) 



-pr c A r X { n ~ 1} - 2 £ 
j'=i 

= -pr c A,£f } - 2 £ 

Z'=l L 



<9 r £n 



(21+1) 



in 



<9 r £n 



<9 r £ri 



a v («-0' 



<9 r £n 



(0 



1 2 



(22) 

Consecutively solving equations beginning from the top, one can obtain S„(r) (21) 
up to an arbitrary order in expansion parameter h/(pr c ). The terms correspond- 
ing to odd I contribute to the amplitude of R n (r), while the terms corresponding 
to even I contribute to the phase of R n (r). 

However, the power of the WKB method reveals fully itself in those cases when 
it suffices to take account for the first two terms, (r) and E^(r), only. It 
should be noted that equation (18), or (20), has its analogue in optics, describing 
the propagation of waves in media [2D]: quotient P n (r)/p corresponds to the index 
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of refraction. Terms E^(r) (I > 2) are negligible, if the variation of the refraction 
index along the ray of propagation becomes appreciable on the distances which 
are much larger than the effective wavelength, 

|p- 1 a-(nP„)|- 1 »n|P n |- 1 , (23) 

where n points in the direction of the ray which in our case is orthogonal to the 
boundary. In view of (19) this yields condition 

« 1, (24) 



\p 2 r 2 — H 2 (n — yu) 2 | 3 / 2 
which is satisfied in the quasiclassical region, 

pr ^ h\n — fj,\, pr ^> h, (25) 
as well as in the deeply nonclassical region 

pr <C h\n — fx\, \n — /i\ >> 1. (26) 
Solving the first two equations in (22), we get 

r 

E?>(r) = ±J drP n (r), E«(r) = -^pr c ln[rP n (r)]. (27) 

It should be noted that the second equation in (22) is nothing more but the 
continuity equation, d • j n , for the partial wave current, 

k = ~ [e^R* n (de in ^R n ) - {de^ W n )e m * R n ] , 

in the quasiclassical region. 

Actually, the WKB approximation (i.e. the neglect of S^(r) with I > 2) is 
efficient in a much more extended regions than those given by (25) and (26). The 
WKB approximation is not valid in the vicinity of r = r t which is defined by 

Pn(r t ) = 0, (28) 

this point is the turning point of a classical trajectory of the particle. In view of 
(27), the two linearly independent solutions to (18) in the quasiclassical region 
are 

RntL(r) = (rP n r 1/2 exp I ±'- J dr P n , r > r t , (29) 



rt 



while the two linearly independent solutions to (18) in the deeply nonclassical 
region are 

R^l{r) = (rn„)- 1 / 2 exp I ±^ / drU n I , r < r t , (30) 




where 

n n (r) = y/h 2 (n-n) 2 r- 2 -p 2 . (31) 

The crucial question is how oscillating solutions (29) in the outer region go over to 
solutions (30) in the inner region and vice verse. To answer this question, one has 

to find a solution in the vicinity of the turning point, where P n (r) = pJ^{r — r t ) 



and II n (r) = p^J jr % (jt ~ T \ and to match it with the solutions in the inner and 

the outer regions. But a more instructive way, as stated in [17], is to consider 
formally R n (r) as a function of complex variable r and to perform transitions 
along paths in the complex plane, where condition (24) is satisfied. Namely, a 
path goes along the real positive semiaxis up to the left vicinity of the turning 
point, circumvents it along a semicircle in the upper or lower half-plane, and then 
goes again along the real positive semiaxis. For the transition from the inner to 
the outer region we have 

n n (r) -)« P n (r)e ±i7r / 2 , (32) 
where the (±) sign corresponds to the path in the upper or lower half-plane, and 

R%L(r) R%L(r)e-^ + i£L(r)e*/* (33) 

-^i _ in( r )' i n contrast to R^^r), is divergent in the inner region. For the transition 
from the outer to the inner region we have 

P n (r) -> n n (r)e ±i ^ 2 , (34) 

and 

" -> (35) 



1 

2 

where the account is taken for the fact that R^^r) is real. 

Boundary condition (16) in terms of partial waves takes form 

[sin(p7r)r<9 r .R n (r) + cos(p7r)i? n (r)]| r=rc = 0. (36) 

In the case r t < r c < r we choose 



Rn(r) 



Ri:Ur) - C n (r c , p)R^Ur) e™'\ (37) 



3 i7r/4 
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where the factor of exp(i7r/4) is inserted for future convenience. Coefficient 
C n (r c , p) is determined from condition (36): 



R ( n,Lt( r c) cot(pvr) + [r d r In R ( n J ut {r 



(-) 



C n (r c ,p) 



R^L(r c ) cot(pvr) + [rd r In <| G J ut (r 



(+) 



|r=r c 



2i 



= exp ^ — — / drP n — 2iarctan 



r c P n (r c )h 1 



cot(p7r) 



1 p^ 



2 P2 (rc ) 



r c > r t - (38) 



note that the explicit form of the turning point (i.e. the solution to (28)) is 
evidently r t = h\n — p\/p. In the case r < r c < r t we choose 



Mr) = R%L(r) - C n (r c ,p) [it£J(r) + cR^ n (r) 
where coefficient C n (r c ,p) is determined from condition (36): 



(39) 



C n (r c ,p) 



(r ) 

1 VinV cj 



<il(^) + ci2 niin(+ )(r c ) 



■x 



cot(p7r) + 


r^lni^r)" 


r=r c 


cot(p7r) + |rd r In 




} 


r=r c 



and constant c is to be fixed before long. Using correspondence rules (33) and 
(35), we continue R n (r) (39) to the case r c < r t < r as 



Rn(r) = R ( n %(r)e-^ + R^U^ 4 - 



-C n (r c ,p) { i [i?S ut (r)e- i7r / 4 - R ( ~Ur)e^ 4 



+ c 



/ 4 J- RC") 



3 i7r/4 



The converging wave, -Ri,out( r )> should disappear from the part with coefficient 
C n (r c , p), cf. (37). This fixes the constant, c = i/2, yielding in the case r c < r t < r 



Rn(r) = R^Ur)e'^ + R^U^ 4 - C n (r c , p)B%L(r)J*<\ (40) 



-iTT/4 , p(-) 



( + ) 



with 



C n {r c ,p) = exp j | y drn„ J x 



x 



cot(pvr) + ± - r c n n (r c )^ 1 i 

+ - exp 



cot(p7r) + 



1 P ^ 

2 H2(r c ) 



+ r c Il n {r- c )h 



-i 2 



drTL 



, r c < r t . (41) 



11 



We can rewrite (37) in the case r t < r c < r as 

R n (r) = R { +Ur)e-^ + R^U^ 4 - [e^ 2 + C n (r c , p)] R^U^ 4 . (42) 

We conclude this section by stating that (40) and (42) give the radial compo- 
nents of partial waves of ip(r, ip) (17) in the quasiclassical region in the WKB 
approximation. 

4 Scattering amplitude and cross section 

Going over to asymptotics prftT 1 — > oo, where 

r 

rP n (r) = pr + 0[h 2 (pr)-\ / drP n = pr - -\n - /x|tt + 0[ft 2 (pr) _1 ], (43) 

rt 

we obtain 
where 



*P(r,ip) =^ \r,cp)+^ c >(r,ip), 



(c)( 



(44) 



and 



p i(prh 1 +tt/4) _^ 

^ (c) (r, <p) = { Yl e in *a n e-^ 1 -^ [e^ 2 + C n (r c , p)] + 

\n-fi\<pr c /h 



pr 



+ e i ^a n e-3l»-^C n (r c ,p)L(46) 

|n— n\>pr c /h ) 

To determine coefficient a n , let us recall the asymptotics of the partial wave 
decomposition of a plane wave, 



3 ipr?i 1 cos tp 



\ g'"y gifarft 1 -ir/4) _|_ e i|n|7r e -i(prn 1 — zr/4) 



A( V ?)e i(pr ^ 1 - 7r/4) + A(y? - 7 r) e -( i ^ 1 --/ 4 )j , (47) 



where A(</?) = (27r) 1 e mip is the angular delta-function, A(ip + 27r) = A (99). 




ne2 



The asymptotics of the plane wave is naturally interpreted as a superposition 
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of two cylindrical waves: the diverging one, e iprh /y/r going in the forward, 
ip = 0, direction and the converging one, e~ iprh / y/r, coming from the backward, 
(p — 7T, direction. Namely the converging cylindrical wave should be present 
without distortions in the asymptotics of wave function tp(r, (p) (44), whereas the 
diverging cylindrical wave is distorted and differs from that in (47) ; actually, this 
is the condition that ip(r, <p) be the scattering state solution. Equating terms 
before e~ iprh / 'y/r in ip(r,(p) to the terms before e~ iprh j y/r in the first line of 
(47), we get 



1 



- v2^ exp [ i y nl ~2 ln ~ i 'T\- (48) 

As a result, we obtain the following expression for the r c -independent part of the 
wave function 



,(o) 



\{prh 1 +7r/4) 



(49) 



where 



ip^ir <p) = e [prh 1 cos ^ e i ^* 3 ~ sgn ^-' 7r ] 



2nh 
pr 



/4) 



(50) 



is the incident wave and 

sin(/i7r) 



fo(p, <P) 



y/ 2-Kpfr 



sin(/i7r) 



i(M" 



^nph- 1 sin(^/2) 



(51) 



is the scattering amplitude 0; the sign function is sgn(w) = ±1 at u ^ 0, \u\ 
denotes the integer part of quantity u (i.e. the integer which is less than or equal 
to u) and it is implied that — ix < <p < n. The squared absolute value of (51) yields 
differential cross section (1). Incident wave (50) differs from the plane wave: the 
distortions are due to the long-range nature of the vector potential outside the 
flux tube. The apparent divergence of amplitude (51) in the forward direction is 
a phantom, because (51) is valid under condition y/pr/h\ sin((^/2)| 3> 1, whereas, 
otherwise, the divergence is absent and the discontinuity in incident wave (50) 
at ip — > ±0 is cancelled, yielding ip(°\r,ip) which is continuous and differentiable 
at ip = 0, see, e.g., (2TJ [22] and references therein. Amplitude fo (51) is of order 

y/r~ c O y/hj (pr c ) and is negligible as compared to the contribution from the r c - 

dependent part of the wave function, which can be presented in the following 
form: 



^ c \r, <p) = [fx(p, V ) + f 2 (p, V ) + f 3 (p, <p)): 



3 i(pr/i _1 +7r/4) 



y/r 



(52) 
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where 

f 1 (p, v ) = —l== J2 e m V(H-l-H)^ (53) 



\n-n\<pr c /h 

1 



^2ixph 

f 2 (p,<p) = -—=L= J2 e m V(H-l-HK a(rc?p)) (54) 

V " \n—/j,\<pr c /h 

/3(p,p) = --_!= e m ^-\ n -^C n (r c ,p). (55) 



|n— n\>pr c /h 



We show in Appendix A (see (A. 22)) that amplitude fa (55), although exceeding 
amplitude /o (51), is still negligible, being of order y/r^O[h 1 ^ 6 (pr c )~ 1 ^]. Ampli- 
tude f 2 (54) is of nonnegligible value which is calculated in Appendix A: 

/2(p, V 5 ) = -^y|sin(v?/2)|exp {— 2ipr c ^~ 1 | sin(v?/2)| + ifj,[(p - sgn(y?)7r] - ivr/4} x 

2pr c h~ 1 \sm 3 ((p/2)\ 



x exp < — 2i arctan 



2cot(p7r)sin 2 ((^/2) - 1 



.(56) 



Amplitude f\ (53) is a finite sum of geometric progression, which is straightfor- 
wardly calculated: 

h(p, <p) = iMsW+fr ™ is "M cosOxtt + s e <p/2) (57) 
y up sm(ip/2) 

in the case 

\pr c h- 1 - [/i] = fc^/T 1 - /x] + M + 1 = s c , (58) 

or 

f lip , „) = iMjn«*» {!l*£±V|M!l cos + (Sc + 1/2W2] _ 

y 7rp [ sm((/?/2J 

-tam>/4) ±i}(59) 

in the case 

[^/T 1 + /i] - [/i] - l - ± i = \pr c h- 1 -n\ + \A + \*\ = Sc. (60) 

Since s c ~ prjir 1 ^> 1, amplitude j\ is strongly peaked in the forward direction 
(actually as a smoothed delta- function), whereas amplitude f 2 is vanishing in the 
forward direction. Therefore, the interference between the amplitudes, f*f 2 + 
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fZfi, is vanishing in the same manner as y/\ip\A(ip) = 0. The differential cross 
section of the scattering process is hence a sum of two terms: 



cboV = l/2(P ' V?)l = 2 



• <P 
sin — 

2 



and 



d<7i 



(l . d = |/i(p,^)r = 4r c A^(v9)cos 2 [(pr cV + - c ^ (2H)- 1 
where we have recalled that \i = e$(27r/ic) _1 and introduced function 

1 sin 2 (yep) 



{-IT <(f< 7r), 



(61) 



(62) 



Any sin 2 '(ip/ '2) 

which at y ^> 1 can be regarded as a regularized (smoothed) delta-function, 
lim AM = A(<p), A y (0) = [ dipA y (ip) = 1 + 0(y- 2 ). 



(63) 



Term (61) which is independent of the self-adjoint extension parameter, the 
magnetic flux and even the particle momentum is obtainable in the framework 
of classical theory, for instance, with the use of the Hamilton- Jacobi equation of 
motion, see Appendix B. Amplitude f 2 (54) can be rewritten as 

f 2 (p, <p) = - ex P H + 2*™(pr c ft- 1 ) - 2u n (pr c H- 1 ,p)} } , 

V r | n _ /j,\<pr c /h 

(64) 

where 

S^^iprjT 1 ) = ^nsgn(n - n)n - ^ ^ dr(P„ - p) - ^pr c = 

oo 

= 7^sgn(n - n)ir - i n {pr c h~ v ) - \ (65) 

is the WKB phase shift, £„(s) and u n (s,p) are given by (A. 2) and (A. 3) in Ap- 
pendix A. The derivative of the WKB phase shift, £$™ KB , can be related to 
the classical reflection angle (equal to the classical incidence angle), 9 Tc — 9^, see 
(B.10) in Appendix B: 



d rWKB | 

8n n J 



n=fj,+a/h 



0r c -0oo + -sgn(a)7r. 



(66) 
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Thus, the condition of stationarity of the phase in (64), 



- - - 2 (^»' KB ) 



(67) 



ra=no 



(£w„ is negligible, see Appendix A), coincides in its form with the classical 
relation between scattering angle ip and reflection angle 9 Tc — 9^, see (B.12) in 
Appendix B, and the classical impact parameter can be introduced as 

b = -h(n - ^p- 1 ; (68) 

note also that the first equation in (22) is equivalent to the equation determining 
the squared r-derivative of the Hamilton- Jacobi action, see (B.5) in Appendix B. 

However, the merits of the WKB method are not restricted to determining 
the WKB phase shift, i.e. to the only description of classical elastic reflection. 
The use of the method has also allowed us to find term (62) which is generi- 
cally independent of the self-adjoint extension parameter but depends on both 
the magnetic flux and the particle momentum. This term describes the purely 
quantum effect of diffraction. It should be noted that the quantum effect is well 
separated in the scattering angle from the classical one: namely, diffraction is in 
the forward direction where classical reflection is absent. Although the range of 
angles where diffraction is to be observed is quite narrow, the whole (i.e. inte- 
grated over — 7r < ip < 7i") contribution of the diffraction effect to the total cross 
section is the same as that of the reflection effect: 

= P = 2r - < 69 > 

az az 

and the total cross section is twice the classical one; the latter is needed for the 
optical theorem to be maintained (see [23| 124"]). 

It should be noted that the AB scattering amplitude (see (51)) and, con- 
sequently, cross section (1) were obtained by imposing the condition of regular 
behavior for the particle wave function at the location of the singular (r c = 0) 
magnetic vortex [2J. However, the regularity condition is not the most general 
one that is required by probability conservation. The most general condition 
in the case when a spatial dimension along the vortex is ignored involves four 
self-adjoint extension parameters (the deficiency index is {2,2}) [251 126]. The 
dependence on the self-adjoint extension parameters enters into terms which are 
added to amplitude f (51). Therefore, the cross section depends on four arbi- 
trary parameters or, if a longitudinal dimension is taken into account, on four 
arbitrary functions of z. If an invariance under rotations around the vortex is 
assumed to be a physical requirement, then the cross section still remains to be 
dependent on two arbitrary functions. The problem of seeking some, if any, sense 
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for this arbitrariness is of purely academic interest, since the case of a vortex with 
r c = has no relation to physics reality. 

Returning to the physical case of a vortex with r c > 0, we would like to em- 
phasize that the classical initial equation (see (B.2) in Appendix B), as well as the 
quantum-mechanical one (see (10) with (13)), involves the vector potential (see 
(12)). However, in classical theory, this potential has no physical consequences, 
as is demonstrated, for example, in Appendix B. One may agree with the authors 
of [12] that "the AB effect was already implicit in the 1926 Schrodinger equation" , 
but it has taken more than eight decades before the physical consequence of the 
vector potential from the Schrodinger equation is at last pointed out, see (62). 

5 Discussion and conclusion 

The AB effect, as well as quantization of the flux trapped in superconductors, 
the Van Hove singularities in the excitation spectra of crystals etc, explicates 
an importance of topological concepts for quantum physics [27]. The AB phase, 
e$/ (he), which is acquired by the charged-particle wave function after the particle 
has encircled a magnetic vortex is certainly invariant as the particle momentum 
varies arbitrarily. Nonetheless, a question always is how to observe this phase in 
real experiments, and what are physical restrictions on the particle momentum in 
these experiments. For instance, interference experiments to verify the AB effect 
involve electrons of energies restricted to the range of 10-100 KeV [8j. Interference 
is immanently a wave-optical phenomenon which is described in quantum the- 
ory indirectly, through reference to wave optics (a quantum-mechanical particle 
possesses de Broglie's wavelength and thus is subject to the laws of wave optics). 
Meanwhile, quantum theory provides a direct description of the AB effect as a 
scattering event. A scattering experiment with particles of the wavelength ex- 
ceeding the vortex radius is hardly realistic, and that is why we are considering 
a scattering event with particles of the wavelength which is less than the vortex 
radius. 

By using the WKB method in the present paper, we obtain the scattering 
differential cross section, da/(dzdip), consisting of two parts: one (61) describes 
the classical effect of elastic reflection from the vortex and another one (62) 
describes the quantum effect of diffraction on the vortex. The effects are equal 
each other, when integrated over the whole range of the scattering angle, see 
(69). It should be noted that, although the amplitude yielding the classical 
effect depends on the parameter of the Robin boundary condition (p), the flux 
of the vortex ($) and the particle momentum (p), see (56), all the dependence 
disappears in the differential cross section, see (61); certainly, this cross section 
can be obtained with the use of purely classical methods (the Hamilton- Jacobi 
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equation of motion). 

The amplitude yielding the quantum effect is generically independent of the 
parameter of the Robin boundary condition, see (57)- (60). This signifies that 
diffraction (to be more precise, the Fraunhofer diffraction, i.e. the diffraction in 
almost parallel rays) is invariant as long as probability is conserved. Owing to 
diffraction the AB effect persists in scattering of real particles, e.g., electrons of 
energies in the range from 10 eV to 100 KeV. The classical effect opens a window 
in the forward direction, where the AB effect is to be observed as a fringe shift in 
the diffraction pattern exhibited by the scattering differential cross section, see 
(62). 

We present the window for the observation of the AB effect on Fig. 
1, where the scattering differential cross section and the scattering angle 
are normalized in such a way that the plot remains actually the same for 
50 < pr c (2-Kh)~ l < oo; recall that the present-day realistic restriction is 
10 4 < •pr c {2 r nK)~ 1 < oo. The classical effect given by (61) is unobservable 
at \cp\ < 107r/i(pr c ) _1 , whereas the diffraction oscillations given by (62) are 

1 1 /3 

unobservable on the classical background at 2 arcsin [10n/(7ipr c )] ' < M < tt; at 
107T^,(pr c ) _1 < \ip\ < 2 arcsin [10^./(7rpr c )] 1 ^ 3 , both effects are unobservable being 
indistinguishable from the |/s| 2 -background which is of order of [h/ '{pr^)] 1 ^ ', see 
(A. 22) in Appendix A. A gate for particles propagating in the strictly forward 
{if = 0) direction is opened when the vortex flux equals an integer times the 
London flux: more than 90% of the diffraction cross section is accumulated in 
a peak centred at <p = and having width 2irh/(pr c ), see the dashed line on 
Fig. 1. The gate for the strictly forward propagation of particles is closed when 
the vortex flux equals a half-integer times the London flux: a dip at (p = 
is surrounded by two symmetric peaks accumulating more than 85% of the 
diffraction cross section (each peak having width 2irh/(pr c )), see the solid line 
on Fig. 1. Thus, the classically forbidden propagation of particles in the strictly 
forward direction is controlled by altering the amount of the vortex flux. 

This scattering AB effect was never observed experimentally, since a detailed 
proposal for its observation was just recently elaborated [28]. Such an experiment 
is in our opinion of fundamental significance and seems to be quite feasible with 
the present-day facilities providing for sources of bright and yet coherent electron 
beams. 
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Appendix A. Calculation of amplitudes / 2 and / 3 

We present amplitude f 2 (54) in the form 

1 

/2(P,¥>) : 



y/2-irph- 1 

x ^2 exp{i[n(f + /2Sgn(n- /2)n-2^ n (pr c h' 1 )-2uj n (pr c h~ 1 ,p)]} ,(A.l) 

\n—^\<pr c /h 

where, see (38), 

Cn(s) = ^ J dr P n = a/s 2 - (n - p) 2 -\n- p\ arccos(|n - /i|/s), (A.2) 

^s 2 - (n - /i) 2 



%-/z|/p 



w„(s, p) = arctan 



(A.3) 



_cot(prr) - L^-(n-^ , 
and we have introduced dimensionless variable s = pr c / ft. To calculate the sum 
in (A.l) in asymptotics s >> 1, we use the Poisson summation formula 

\n— /i|<s 



ir;(n,s) 



= J2 [ dnexp{i{ri(n,s) - 2iml]} + \e^ s +^ + )-e^- s ~> s \ (A.4) 
,_™ J z z 



where s± = [s ± //]. If s + + s_ ^> 1 and 77(71, s) is convex upwards, 8 l^l'^ < 0, 



dn 2 



on interval — s_ < n < s+, then only a finite number of terms in the series on the 
right-hand side of (A.4) contributes to the leading asymptotics at s >> 1, and one 
can use the method of stationary phase for its evaluation. Namely, if equation 

d 



dn 



rj(n, s) 



2ttIj = (A. 5) 



determines a stationary point inside the interval, — s_ < rij < s + , for some values 
of / denoted by lj, then 

1/2 



= ]Texp {ifofa, s) - 2™^]} 

\n—/j,\<s lj 



2ne~' lnl2 



+ 0(1). 
(A.6) 
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In the present case we have 

r)( n , s) = nip + /isgn(n - p)n - 2[f n (s) +w n (s,p)]. 
Differentiating (A. 2) and (A. 3) over n, we obtain 
d 



^^n(s) = -sgn(n - p) arccos(|n - p\/s) 



(A.7) 



(A.8) 



and 



9 ( \ 





1/2 






cot(p7r) - § j 


{[l-f 2 ?)" 


cot(p7r) - \ > + s 2 


;i-(n?) 2 f 



(A.9) 



(A. 9) is negligible as compared to (A.8), being of order 0(\n—p\/s 3 ) at |n— /i| <C s, 
s>l, where the WKB approximation is valid, see (25). Therefore equation (A. 5) 
takes form 

ip + 2sgn(nj — p) arccos(|nj — p\/s) — 2nlj = 0. (A. 10) 

Choosing the range for p as — 7r < p < it, we find that a solution to (A. 10) exists 
at lj = only. Denoting the solution by rij = n , we obtain 



and, hence, 



n - p = -s sgn(y?) cos(</?/2), 
<9 2 , 



<9n 2 



r/(n,sj 



s| sin(</?/2)| 



(A.ll) 
(A.12) 



We see that the term corresponding to / = in the right-hand side of (A. 4) is of 
order 0(^s), whereas all other terms are of order 0(1). As a result, we obtain 
expression (56). 

Amplitude f 3 (55) is presented in the form 



fs(p,p) 



^2nph 



[n</3+/^sgn(?i-/i)7r] 



X 



\n—fi\>s 



x exp 



2\/(n — p) 2 — s 2 — 2|n — yu|arccosh(|n — p\/s) 



x 



x 



cot(p7r) + \ 



2 (n-ti) 2 -s 2 



\/(n- p) 2 - s 2 



cot(pTr) + | (n _; } 2 2 _ s2 + v/(n-p) 2 -s 2 



+ -exp 



2a/ (n — p) 2 — s 2 — 2(n — /i)arccosh ^ 



|n — yU| 



(A.13) 
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where (41) is taken into account. The sum in (A. 13) at ip = is evaluated 
asymptotically at s ^> 1 by converting it into integral: 



2r c 



71 



s COS 



(mtt) [lW(s)-ilW(s)]+0(l/Vt)}, (A.14) 



where 



I (s) = 4 / dvexp 2s (farccoshu — \/v 2 — 1 



x < 1 + 4 exp 



4s (i? arccoshf — \/v 2 — 1 



cot(pvr) + i^23i - sVf 2 - 1 



cot(/wr) + |^2^j + s\/v 2 - 1 
cot(pvr) + |^ - svV - 1 



cot(pTr) + + sVv 2 - 1 



(A.15) 



and 



x < 1 + 4 exp 



4s (v arccoshf — y/v 2 — 1 



CO 

/(2)( s ) =2 /" dt; 



cot(p7r) 



11 

2? 



cot(pvr) + + sy/v 2 - 1 



-i 



(A.16) 



In the cases of |cot(/?7r)| 3> s and |cot(/97r)| ~ s, the integrand in (A.15) is 
maximal at v = 1, and one can use the Laplace method for the estimation of 
/W(s) in this case, see, e.g., [29]. Namely, let g{v) be continuous function with 
< g(l) < oo and h(v) be continuous differentiate function obeying conditions 
£h(v) < at v > 1 and « -a(v - l)^ 1 at v ->■ 1 (v > 0), then 



i> V f / Vas- 



(A.17) 



where T(y) is the Euler gamma-function. In the present case g(l) = 1, h(l) = 0, 
a = 2\/2, v = 3/2, and we obtain 



(A.18) 



A similar estimate is obtained for I^ 2 \s) (A.16) in the case of | cot(p7r)| ^> s: 



/(2)(S) = I r (I) ^ 



2\-l/3 



(A.19) 
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In the case of |cot(p7r)| <C s, the integrand in (A. 16) is maximal at v = vq 
where vq = 1 + 2~ 5 / 3 s~ 2 / 3 . The integrand in the vicinity of the maximum is 
approximated as 

2exp {-9(2s) 4/3 exp[10(2s)~ 2/3 ](w - v ) 2 } , 
and the integral is estimated as 



(A.20) 



The integrand in (A. 15) in the case of | cot(/wr)| <C s changes sign at v = v . 
Similarly to (A.20), we obtain the estimate for integral I^'(s) in this case: 



7 (D / 



1 / 1 \ 2/3 
~3 \2^) 



(A.21) 



In the case of | cot(p7r)| ~ s, the numerical analysis of I^ 2 '(s) yields an estimate 
of order 0(s~ 2 / 3 ). Hence, we obtain the estimate for amplitude f 3 (p,Lp) (A. 13): 



where constant a is independent of s, [i and (p. 



(A.22) 



Appendix B. Hamilton- Jacobi equation of motion 
and classical cross section 

In classical theory, we consider scattering of a point charged particle by an impen- 
etrable tube containing magnetic flux $. The Hamilton function corresponding 
to the motion outside of the tube is 



H 



2m 



P T + -oiPe- ft) +P Z 



dS 1 



where // = e$(27rc) l , and the Hamilton- Jacobi equation (see, e.g., 

dS\ 2 }_(d^_ A 2 /^Yl-n 

A general solution to (B.2) is 

r 6 z 

S = -Ht + / drp r + / dOp e + / dzp z . 



(B.l) 



is 



(B.2) 



(B.3) 
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Primary conservation laws are 

1 

p e = a, p z = a, — 

2m 



2 / A 2 2 

Pr + ^(Pe-V) +P Z 



E (B.4) 



with a, a and E being arbitrary constants (E > 0). Taking the conservation laws 
into account, we obtain action S (B.3) in the form 

r 

S = -Et + J dr v / 2mE -a 2 - r~ 2 (a - /i') 2 + a0 + az + const. (B.5) 

Secondary conservation laws are 

dS n dS n dS n ,„ „, 

with /3 Q , /3 5 and /3e being arbitrary constants. Thus we obtain relations 

fia = e-f dr - CtzJtH (B.7) 
J y/2mE -a 2 - r~ 2 (a - n') 

r 

^ ~ Z ~ I ^ ^2mE -a 2 - r-\a - fi') 2 ' ^ 

r 

/771 
dr (B.9) 
^j2mE — a 2 — r 2 (a — fi') 2 

which determine a trajectory of the particle, i.e. the functional dependence of 9, 
z and t on r. One can conclude that the classical trajectory is independent of the 
enclosed magnetic flux, since // can be absorbed into constant a: a — // — » <x 
The trajectory is symmetric with respect to the point of reflection, r = r c . 

As a consequence of (B.7), we obtain the following expression for the incidence 
angle which is equal to the reflection angle: 

TV 

r~ 2 a ( a 



rc -0oo= dr - = = - arcsin — , (B.10) 

' y/p 2 - r- 2 a 2 \pr c ' 



where p = \j2mE — a 2 {2mE > a 2 ). Defining the impact parameter as 

b = -ap^ 1 (—r c <b <r c ) (B.H) 

and using elementary tools of plane geometry, we determine the angle of deflection 
(called also the scattering angle) in range — n < <p < n as 

p = sgn(6)7T - 2(6 r . c - 6^) = 2sgn(6) arccos(|6|/r c ). (B.12) 

Hence, the scattering differential cross section is 

^ class db r c . LD. . 

-|—|— = -^~ = 77 sin £ . B.13 

dzdp dip 2 1 2 1 y J 
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